Quantum algorithms for spin models and simulable gate sets for quantum computation 
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We present elementary mappings between classical lattice models and quantum circuits. These 
mappings provide a general framework to obtain efficiently simulable quantum gate sets from exactly 
solvable classical models. For example, we recover and generalize the simulability of Valiant's match- 
gates by invoking the solvability of the free-fermion eight-vertex model. Our mappings furthermore 
provide a systematic formalism to obtain simple quantum algorithms to approximate partition 
functions of lattice models in certain complex-parameter regimes. For example, we present an 
efficient quantum algorithm for the six-vertex model as well as a 2D Ising-type model. We finally 
show that simulating our quantum algorithms on a classical computer is as hard as simulating 
universal quantum computation (i.e. BQP-complete). 



1. Introduction. — Interrelations between quantum 
information theory and statistical mechanics have re- 
cently proven to be particularly fruitful. Indeed, it turns 
out that several connections between these fields exist, 
such that techniques from either one of these research 
areas may be used to gain insights in the other [H, 0, S 0] ■ 

A natural context to study such connections to sta- 
tistical mechanics regards the investigation of the com- 
putational power of quantum computers — in particular, 
understanding their potential and limits, and the bor- 
derline between classical and quantum computation. Re- 
cent findings indeed point to a fruitful cross-fertilization 
between this subject and statistical mechanics. For ex- 
ample, in f2] it was proved that the exact solvability of 
the 2D Ising model can be invoked to prove that the 
"toric code states" @, regarded as resources for quan- 
tum computation, do not yield any speed-up with respect 
to classical computation. This indicates that, under cer- 
tain circumstances, exactly solved models may be related 
to the efficient classical simulation of quantum computa- 
tion. In other work Q , efficient quantum algorithms 
are obtained to approximate difficult quantities such as 
the Potts model partition function frj . This suggests that 
quantum computers may become useful tools to study 
problems in statistical mechanics. 

Despite the existence of these interesting results, such 
findings typically occur as isolated mappings, and a sys- 
tematic approach seems to be missing. Indeed, although 
recent algorithms such as 0] are significant achievements, 
they require involved arguments and are not easily gener- 
alized. One may therefore hope that quantum algorithms 
can be constructed to tackle statistical mechanics prob- 
lems in a simpler and more systematic way. Furthermore, 
in statistical mechanics a lot of effort has gone in develop- 
ing techniques which aim at solving spin systems, such 
as the Jordan- Wigner transform and the Yang-Baxter 
equation Q. It would be highly beneficial if clearcut 
mappings existed which would allow one to use such es- 



tablished techniques to study the computational power 
of quantum computers in a general setting, going beyond 
results such as e.g. 

We believe that in the present paper significant steps 
are taken towards achieving the above goals: 

(i) We present elementary, general constructions that 
relate arbitrary classical lattice models with quan- 
tum circuits. 

(ii) Using these mappings, we show that from any clas- 
sical model where the partition function can be ef- 
ficiently evaluated, we obtain a class of quantum 
circuits that can be efficiently simulated classically. 

(iii) These mappings also immediately provide a simple 
framework to construct quantum algorithms that 
approximate the partition function of these models 
in certain — typically complex-valued — parameter 
regimes. We furthermore show that approximating 
the partition function of certain models is equiv- 
alent to simulating quantum computation alto- 
gether, i.e., these models are BQF-conrplete [11]. 

More concretely, these general findings allow us to, e.g., 
recover and generalize a construction of Valiant [9^ , yield- 
ing classes of simulable gate sets from the solvability of 
the 8- vertex and 32- vertex model with free-fermion condi- 
tion. We furthermore obtain a class of simulable gate sets 
by invoking the simulability of the 2D Ising model with- 
out external field. We will also show that the 6-vertex 
model as well as a 2D Ising-type model are, in certain 
complex parameter regimes, BQF-complete models. 

2. Quantum circuits as vertex models. — Next 
we show how a general class of classical spin systems, 
called "vertex models" , is naturally associated with quan- 
tum circuits. A vertex model (VM) consists of classical q- 
state spins Se G {0, 1, . . . , q~ 1} placed on a lattice, where 
a spin Se is associated with each edge e of the lattice, and 



2 



where (many-body) interactions occur in the vertices of 
the lattice. The energy of a spin configuration s = (se) 
is given by a sum of local contributions: H{s) — J^a 
Here h"" is the interaction at site a, which is a function 
of those spins situated at all edges incident with a. For 
example, in the case of a VM on a 2D square lattice, 
h" = h''{su, Sd, si, Sr) is a 4-body interaction depending 
on the spins "up" , "down" , "left" and "right" of vertex 
a (see Fig. [l]). 

Instead of defining an energy function H, a VM is 
frequently expressed directly in terms of its Boltzmann 
weights W{s) — e^'^^^*) (where /? denotes the inverse 
temperature). This typically amounts to prescribing 
which spin configurations are allowed to have a nonzero 
weight. Note that W{s) = JJ^ w"(s) is a product of local 
weights w'*(s) = e~^'^ For example, the "eight- vertex 
model" is a model on a 2D square lattice with two-state 
spins, where only eight out of the 2'* = 16 spin configu- 
rations around each vertex are allowed to have a nonzero 
weight. 

A central quantity for any spin system is the partition 
function 

Z = ^e-^«(«), (1) 

since all thermodynamical properties of the system can 
be derived from Z. 

We now show how vertex models can be mapped to the 
quantum circuit model. We will illustrate the mapping 
for a q-state vertex model on a 2D tilted square lattice 
(see Fig. [T]). The construction is, however, not limited to 
this case, as we will see below. We choose open bound- 
ary conditions and furthermore fix the spins at the left 
and right boundary of the lattice (where one has, say, N 
spins on each side) in some arbitrary configurations de- 
noted by i = (sf , . . . , s^) and R = (sf , . . . , sfj). In the 
horizontal direction, the lattice contains M = poly(7V) 
sites. Note that, as every vertex of the lattice has degree 
4, there are possible spin configurations at each site, 
and therefore possible Boltzmann weights w"'{i,j, k, I); 
here i,j,k,l denote the four spins at the site a. To 
simplify notation we will denote the Boltzmann weights 
by w]^i, hence omitting the dependence on a (but keep- 
ing in mind that Boltzmann weights are generally site- 
dependent). We will organize these Boltzmann weights 
in a X matrix denoted also by , by grouping the 
spins i,j {k,l) to be the row (column) index of the ma- 
trix: w"" = ^w^^i\ij){kl\. For example, in the case of the 
8VM P, w° is a 4 X 4 matrix of the form 
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The six-vertex model [8| is obtained by additionally set- 
ting u;00 = = «;ii. 




Figure 1; [Top left] Vertex model on a tilted 2D square lattice. 
On every edge, a g-state classical spin (illustrated in red) is 
situated. At every vertex, an interaction takes place between 
the four incident edges i,j,k,l, with corresponding Boltzmann 
weight The spins on the left and right boundary edges 

are fixed in certain configurations L and R (indicated in blue) . 
[Top right] Corresponding quantum circuit C, acting on a lin- 
ear chain of quantum systems. The g* Boltzmann weights 7j;^-'j 
at each site constitute the matrix elements of a quantum gate 
acting on two neighboring systems (represented by a box). 
The boundary conditions L and R are mapped to basis states 
jL) and \R). [Bottom left and right] Vertex model on triangu- 
lar lattice and corresponding quantum circuit built of 3-local 
gates. 

The connection to the quantum circuit model is imme- 
diately obtained by considering the matrices as 2-local 
(possibly non-unitary) quantum gates acting on pairs of 
nearest-neighbor g-level quantum systems. More partic- 
ularly, we consider a linear chain of g-level systems, 
and a quantum circuit C — Yia acting on this system 
as depicted in Fig. [TJ It is then an easy exercise to show 
the following identity: 

Z^i^={L\C\R). (3) 

Here Z^^^ is the partition function of the VM with 
boundary conditions L and R, and \L) and \R) are A^- 
party computational basis states associated with the spin 
configurations L and R in the obvious way. 

Equation ^ establishes in a simple way a direct rela- 
tion between VMs and the quantum circuit model. The 
mapping is akin to well known mappings from 2D classi- 
cal to (H-l)D quantum systems; the fact that we are 
considering a tilted 2D lattice, allows one to read off 
the quantum circuit C immediately. A number of con- 
clusions can immediately be drawn from this correspon- 
dence. This is shown in the next two sections. 

2.1. Simulable gate sets from solvable VMs. — 

First, it is clear from Eq. ^ that every VM for which 
the partition function Z^^ can be evaluated (or approx- 
imated) efficiently, gives rise to a quantum circuit for 
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which the matrix element {L\C\R) can be evaluated ef- 
ficiently classically. We will illustrate the power of this 
mapping by giving three non-trivial examples of simu- 
lable gate sets which can be obtained in this way. The 
first example is obtained by considering the 8VM. It is 
known that the partition function Z^^\^ of this model, 
with possibly complex Boltzmann weights, can be evalu- 
ated efficiently for every L, R and for every system size 
N , as long as the Boltzmann weights ([2]) additionally 
satisfy the condition 
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Under this condition, the model can be mapped to a 
system of non-interacting Fermions [To| . which can be 
solved analytically. The correspondence ([3]) immedi- 
ately implies that all circuits C built out of gates of 
the form with the additional condition ([4]), can be 
simulated efficiently classically in the sense that {L\C\B) 
can be evaluated efhciently for all basis states \L) and 
\R) . Interestingly, this recovers a well-known result Q, 
where it was also showed that the above gates (which 
are called "match-gates") are efficiently simulable clas- 
sically. The present considerations show that the sim- 
ulability of match-gates is implicit in the solvability of 
the 8VM model, a result which has been known since 
long. Also the connection between match-gates and free- 
fermion systems, which was (re)discovered in the context 
of quantum computation in ^13| | , has been established for 
some time in the context of the 8VM (see also [Tl|). 

As a second example, we consider a VM on a 2D tilted 
triangular lattice as in Fig. [U with boundary conditions 
L and R. Going through a similar argument as above, one 
finds that the partition function of this model can 

again be written in terms of a quantum circuit C = lu" 
as in ([3]), but where now the gates w° act on triples rather 
than pairs of g-level systems. Indeed, since every site 
in the lattice has degree 6, the Boltzmann weights are 
organized in a g'^ x q'^ matrix, such that each elementary 
quantum gate is a 3-local gate. The structure of the 
circuit is also given in Fig. [1] We now recall that a 
particular q = 2 VM of this form, the "32-vertex model" 
is known to be solvable [l2i | (also by using a mapping to 
free fermions). Because of space limitations we omit the 
specific form of the 32VM Boltzmann weights, and we 
refer to e.g. [T^ . 

As a third example, we remark that there exist q — 2 
solvable (free-fermion) VMs on a 3D square lattice, see 
e.g. [l^. A 3D model on a tilted lattice would correspond 
to a simulable quantum circuit acting on a 2D array of 
qubits. Note that, as the degree of the 3D lattice is 6, 
the elementary gates in the circuit are 3-qubit gates. 



can be efficiently approximated on a quantum computer 
using the "Hadamard test" (see e.g. @). Using ([3]), 
this immediately yields efficient quantum algorithms to 
approximate partition functions of VMs in those param- 
eter regimes for which the matrices are unitary. The 
type of approximation which can be achieved by the 
Hadamard test, is an "additive approximation". This 
entails that this quantum algorithm returns as a result a 
number c in poly(Af) time which satisfies 

\c~{L\C\R)\<piNr\ (5) 

with a success probability that is exponentially (in N) 
close to 1, where p{N) is some polynomial. We note that 
this type of approximation is achieved in many recent 
quantum algorithms [1, [l^. As an example of the 
present technique, an efficient quantum algorithm to ap- 
proximate the partition function of the 8VM (and, a for- 
tiori, the 6VM) is obtained for all Boltzmann weights for 
which the matrices ^ are unitary. Note that no restric- 
tion is required regarding translation invariance of the 
model: the w°'^s may be different at each site, as long as 
these matrices are unitary. On the other hand, it is clear 
that, in order for w° to be unitary, some of the Boltzmann 
weights have to be negative (or even complex) numbers. 
Therefore, it seems difficult to obtain with this technique 
approximations of partition functions with positive and 
thus "physical" Boltzmann weights. This problem seems 
to be present also in other work, see e.g. [3|, IJ], and it is 
at present unclear whether there exist a general method 
to overcome it. 

We now show that our quantum algorithms are indeed 
non-trivial, by considering the computational complexity 
of the problems that can be solved with these algorithms. 
In particular, we will show that approximating the parti- 
tion function of the 6VM in the abovementioned unitary 
parameter regime is as difficult as simulating all quan- 
tum computations, and hence unlikely to be possible in 
an efficient way classically. More precisely, we will show 
that approximating the partition function Zg^^^ of the 
6VM on a (tilted) 2D square lattice, with staggered left 
and right boundary conditions L — R — (0101 . . . ), and 
where the matrices are unitary, is a BQP-complete 
problem [l!]. To achieve this, we need to show that every 
quantum computation can be reduced, with polynomial 
computational effort, to the approximation of Z^^^ as in 
([5]). To prove this, we use that the exchange interaction 

^^ox ^ CTx ® CTx + (7y ® (7y + (7^ 1^ (7z (6) 

is (encoded) universal for quantum computation p^ . 
The corresponding unitary gates U — e^*^°'^ are of the 
form ^ with non-zero entries 



2.2. Quantum algorithms for VMs. — We now 

consider a second application of identity ([3]) , namely the 
construction of quantum algorithms to approximate par- 
tition functions of VMs. Also here the argument is par- 
ticularly simple. The key point is that, given any poly- 
sized unitary quantum circuit C, the quantity {L\C\R) 
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WqI = ism[2t). (7) 

The required encoding uses four-qubit logical states of 
the form |0l) cx (|01) — |10))®^. From this encoded uni- 
versality, it can be shown (see e.g. 6]) that the problem 
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of providing an additive approximation (as in ([5])) of the 
number (Ol|®^C|Ol)®^ is BQP-complete, where C is an 
arbitrary poly-sized quantum circuit acting on A'^ logical 
qubits, and where each elementary gate in C has the form 
gitHox^ In order to connect this problem to the 6VM par- 
tition function, we first note that the unitary two-qubit 



gates e 



itHn 



in C correspond to Boltzmann weights w'^ of 



a 6VM. Second, we need to be careful with the boundary 
conditions: since |0l) is an entangled state, we cannot 
simply interpret it as a boundary condition as in ([3]), 
since \L) and \R) need to be computational basis states. 
However, we note that |01) — |10) can be prepared from 
a basis state |01) with a unitary gate V which "obeys" 
the 6VM structure: simply take V to have the structure 
with non-zero elements 
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This allows one to write 



{OLr'c\OLr" = {oi\^'''c'\oiY^''\ (9) 

where C = [y^(S2Nqy(S2N ^ -y^g ^^^^ identify 
(01|®2WC'|01)®2^^ with the partition function of a 6VM 
with staggered left and right boundary conditions. This 
proves the BQP-completeness of the 6VM. 
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Figure 2: [Left] Edge model on a 2D square lattice. On ev- 
ery vertex, a g-state classical spin (illustrated in red) is sit- 
uated. Along every edge, an interaction takes place between 
the two spins i and j at the incident vertices, with corre- 
sponding Boltzmann weight Wij. The spins on the left and 
right boundary vertices are fixed in certain configurations L 
and R (indicated in blue). [Right] Corresponding quantum 
circuit C, acting on a linear chain of quantum systems. The 
Boltzmann weights Wij at each horizontal edge constitute 
the matrix elements of a quantum gate acting on one q-level 
system; the weights at each vertical edge constitute the 
diagonal entries of a diagonal two-local gate acting on two 
neighboring systems. The boundary conditions L and R are 
mapped to basis states \L) and \R). 



3. Edge models and quantum circuits. — In anal- 
ogy with the above results for VMs, next we present map- 
pings between "edge models" and quantum circuits. In 
contrast with a VM, an edge model (EM) is a spin sys- 
tem on a lattice where g-state spins Sa reside on the sites 
a of the lattice, and two-body interactions h'^{sa, Sb) oc- 
cur along the edges e = a6 of the lattice, between pairs 
of spins. Important examples of EMs are the Ising and 
Potts models 0- We will illustrate our results for an 
Ising-type model with 2-state spins arranged on a 2D 
rectangular lattice (see Fig. but the results will be 
easily generalized. The energy function of such a model 
is H{s) — J2e ^'^i^a, Sb) (where the sum is over all edges 
e = ab) and we denote the corresponding Boltzmann 
weights by wf^- = exp[— /3/i^(i, j)] (with i,j = 0,1). No- 
tice that, by allowing complex coupling strengths, the 
Boltmann weights wfj can take arbitrary (complex) val- 
ues. The dimensions of the lattice are N x M, where 
M = poly(7V) is the horizontal dimension, as before, but 
this time we do not consider a tilted lattice. Finally, we 
fix left and right boundary conditions L and R. 

We now associate with this spin system a quantum cir- 
cuit in the following way. Consider a linear chain of N 
qubits and a quantum circuit C consisting of AI "layers" 
of one- and two-qubit gates, which has a structure as de- 
picted in Fig. [21 There are two types of layers in C, which 
alternate each other. The first type of layer is associated 
with a slice of horizontal edges in the lattice; it consists of 
N single-qubit gates acting on the N qubit wires, where 
each gate has the form J2i j=o ^ljl*)OI- That is, the four 
possible Boltzmann weights along an edge form the en- 
tries of a single-qubit gate. The second type of layers is 



associated with a vertical slice of edges, and consists of 
A'^— 1 nearest- neighbor, diagonal (and thus mutually com- 
muting) two-qubit gates of the form j=o ^Ij I*-?) 
Thus, here the Boltzmann weights are the diagonal ele- 
ments of these diagonal gates. With this circuit C, one 
finds that 

z4^ = (i|c|i?). (10) 

Here is the partition function of the EM with 

boundary conditions L and R, and \L) and \R) are de- 
fined as before. Thus, by suitably identifying the Boltz- 
mann weights of a model with single- and two-qubit uni- 
tary operations, the partition function of any EM is as- 
sociated with a quantum circuit. 

3.1. Quantum algorithms for edge models. — As 

in the case of VMs, Eq. PH)) now immediately leads to 
a simple quantum algorithm, via the Hadamard test, to 
approximate in those parameter regimes for which 

the single-qubit gates '^wlj\i){j\ and two-qubit gates 
'^'Wij\ij){ii\ are unitary. Note that a unitary two-qubit 
gate is obtained if and only if \wlj \ = 1. 

Moreover, it is again possible to prove that these al- 
gorithms are associated with BQP-complete problems. 
To show this, one remarks that universal quantum com- 
putation can be achieved by quantum circuits of the 
form considered here. In fact, universality is already ob- 
tained when considering a suitable discrete set of unitary 
gates — corresponding to a discrete set of (complex) cou- 
plings in the corresponding EM. Indeed, it is well-known 



5 



that the single-qubit gates 

h = |0>(0| + |1>(1| (identity), 

H (X (Hadamard), 

P = |0)(0| +e"/«|l)(l| (phase gate), (11) 

together with the two-qubit identity and nearest- 
neighbor controUed-phase gate 

CP = ^(-l)^-^|zj)(zj|, (12) 

are sufficient to achieve universal quantum computation. 
This imphes that the approximation (in the sense of ([5])) 
of an Ising-type partition function on a 2D square 

lattice, with possible Boltzmann weight associated with 
the above universal gate set, and with boundary condi- 
tions L = R = (000 . . . ) (where all boundary spins are 
pointing in the same direction) , is a BQP-complete prob- 
lem. 

Note that the BQP-completeness of this 2D Ising-type 
model (and hence the existence of an efficient quantum 
algorithm to approximate the problem) is remarkable, as 
it is known that the exact evaluation of the 2D Ising spin 
glass in the physical parameter regime is associated with 
NP-complete problems (T7| . 

We finally note that a similar quantum algorithm as 
above can be obtained for the partition function of the 
standard 2D Ising model with magnetic fields, although 
in this case one needs to take a number of pre-factors 
into account in the relation (jTU]) between the partition 
function and the quantum circuit. 

3.2. Simulable gates sets for edge models. — 

By considering exactly solvable EMs, such as the 2D 
Ising model without magnetic field, the mapping (jlOp 
allows one to directly read off simulable gate sets. For 
the 2D Ising model without external field, which has 
an energy function satisfying h'^{Q,0) = h^{l,l) and 
h'^{0,l) = /i'^(l,0), the corresponding single-qubit gates 
of the form expiiaax) , and the two-qubit diagonal gates 
of the form exp(i/3(Tz (g) az)- Invoking that the partition 
function of such a 2D Ising model can be efficiently eval- 
uated (the latter can again be proved by mapping this 
classical system to a quantum system of non-interacting 
fermions), we can conclude that every circuit C built out 
of (Ta;-rotations and nearest-neighbor cr^ (8) cr^- rot at ions, 
can be efficiently simulated classically (in the sense that 
{L\C\R) can be evaluated efficiently). 

4. Summary and outlook. — We have established 
simple but general mappings between quantum circuits 
and classical vertex and edge models. This framework 
has been used (i) to derive simulable gate sets for quan- 
tum computation from exactly solvable classical models; 
and (ii) to derive quantum algorithms to approximate the 



partition function of classical models in certain param- 
eter regimes, showing that approximating certain parti- 
tion functions is as hard as simulating universal quantum 
computation. Our results are not limited to the examples 
given in this paper, and can easily be extended. 

The investigation initiated here is by no means com- 
plete, and faces several interesting challenges. 

Regarding (i), we note that all solvable models pre- 
sented here are solvable via mappings to free fermions. 
The major challenge is to also incorporate in our frame- 
work the powerful formalism of the Yang-Baxter equation 
(YBE) and the Bethe ansatz in order to study the clas- 
sical simulation of quantum computation. In this con- 
text one needs to be careful to take into account that 
the YBE typically provides solutions for classical mod- 
els in the thermodynamic limit, and in a limited regime 
of parameters, and often for systems with a high degree 
of translational symmetry. These limitations imply that 
the YBE cannot trivially be incorporated in the present 
framework, and some additional work is required; this is 
an ongoing investigation. In this context, it is also impor- 
tant to remark that in this paper we have only considered 
classical models for which the partition function can be 
evaluated exactly. However, for an efficient simulation 
of the corresponding quantum circuits, an approximate 
solution would be sufficient. Our current investigations 
indicate that such approximations may play a role when 
employing the YBE in this study. 

We also mention that many of the simulability results 
in this paper are not restricted to unitary quantum cir- 
cuits: efficient classical simulation of certain non-unitary 
circuits can be achieved in a similar way (describing e.g. 
post-selected quantum computation). 

Regarding (ii), we believe that the main merit of the 
obtained quantum algorithms is twofold. First, prov- 
ing BQP-completeness of a spin model, as we have done 
for the 2D six-vertex model as well as a 2D Ising-type 
model, strongly suggests that "solving" such a model is 
a computationally difficult task. We may therefore learn 
something about the complexity of spin systems via the 
present results. Second, BQP-complete problems by def- 
inition capture what the power of quantum computation 
is. By relating quantum computation to well studied 
classical spin systems, we believe that new insights may 
be gained about the power of quantum computers. 

The main challenge in the context of these quantum 
algorithms is clear: it would be of great interest to ob- 
tain quantum algorithms which can treat classical spin 
systems in their "physical" parameter regime, i.e., con- 
sidering real and positive Boltzmann weights. This would 
make it possible to use a quantum computer to simulate 
classical spin systems. We are currently studying this 
issue within the framework presented in this paper. 
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After completing this work, we became aware of other 



recent work regarding quantum algorithms to approxi- their investigation is related to our work 
mate partition functions of classical spin systems p^ : 
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